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I. The reactive energy density TZ(r, t) 

We give a physically compelling definition of the instanta- 
neous reactive energy density associated with an arbitrary time- 
domain electromagnetic field in vacuum [1]. In Heaviside- 
Lorentz units, where sq = /io = 1, it is given in terms of 
the energy density U(r,t) and the Poynting vector S(r,t) by 

K(r, t) = CD 

This is a field- theoretic version of the rest energy of a rela- 
tivistic point particle with total energy E and momentum p, 

E = y / E 2 -C 2 p 2 . 

We may interpret (1) as follows: at space-time points (r,i) 
where |S| < U, the energy flow is insufficient to carry away 
all of the energy in the form of radiation. The (momentarily) 
abandoned 'rest' energy is reactive. 

In terms of the electric and magnetic fields (E,B), we have 

U =\(E 2 + B 2 ), S = ExB 
and 1Z reduces to the simple expression 

1Z = yj '\(E 2 - B 2 ) 2 + (E • B) 2 > 0. (2) 

This shows that at each space-time point (r,t) we have 

1Z = & E 2 - B 2 = and E • B = 0, (3) 

which are precisely the conditions for a pure radiation field. For 
a generic EM field, 1Z is strictly positive almost everywhere 1 
in space-time and approaches zero, as it must, only in the far 
zone. Fields for which 1Z vanishes identically, called null 
fields, consist of pure radiation. The simplest null fields are 
traveling plane waves. An interesting example of null fields 
with sources, resembling a spinning black hole in general 
relativity, was constructed in [2]. It was this example that 
inspired the general study of reactive energy density in [1]. 

Just as the rest energy Eq defines the mass m of the point 
particle by Eq = mc 2 , so does 1Z define the electromagnetic 
inertia density X by 

K(r,i) = l(r,t)c 2 . 
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^ere almost everywhere means that 7Z can vanish only on lower- 
dimensional hypersurfaces of space-time. If E 2 — B 2 and E • B are 
independent, (3) implies that 1Z = on a 2D space-time surface whose 
time slices are, in general, time-dependent curves in space. For the standing 
plane wave in Example 2 below, E • B = 0, so (3) reduces to one condition 
and 1Z vanishes on the traveling planes (6), which form 3D hypersurfaces in 
space-time whose time slices are snapshots of the planes at a given time t. 



Whereas m and Eq measure impedance to acceleration, X and 
1Z measure impedance to radiation. Like E , 7Z is Lorentz 
invariant, i.e., it has identical values in all uniformly moving 
(inertial) coordinate frames. For narrowband fields, the time 
average of 1Z is expected to reduce to the known, stationary 
reactive energy density. Thus 1Z is a transient or 'ultra- 
wideband' version of the latter, local in time as well as space. 

We compute 1Z(r, t) explicitly for two fields representing the 
extremes of space-time localization: 

1) A general time-dependent electric dipole field. This is 
local in space-time. 

2) A standing plane wave obtained by adding two plane 
waves of frequency uj > traveling along ±z. This is 
localized at two points in the 4D frequency-wavenumber 
domain, hence highly nonlocal in space-time. 

In Example 1, we find that the reactive energy oscillates 
around the dipole, as shown in Figure 1, and decays to zero 
in the far zone as expected. 




Fig. 1. For a frequency-modulated Gaussian electric dipole, the near- 
field pattern of the reactive energy density oscillates between the two 
forms shown above; see [1] for details. 

In Example 2, we have 

U = E 2 [cos 2 (kz - ujt) + cos 2 (kz + ut)] 
S = zE 2 [cos 2 (kz - ojt) - cos 2 (kz + ujt)}, 

where k = uo/c and E is the amplitude of the electric fields 
of the traveling plane waves. This gives 

1Z = 2E 2 1 cos(kz - ujt) cos(kz + ut) \ . (5) 



Thus 1Z vanishes on the traveling nodal planes z = zf(t\ 
where 



zf(t) = ^l±l)l ± ct e = o, ±1, ±2 
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(6) 



and 7£ > elsewhere. This is shown in Figure 2. 




Fig. 2. The reactive energy density lZ(z,t) of the standing plane 
wave in Example 2, showing the moving nodes (6) where 1Z — 0. 



The two plane waves traveling along ±z are null, i.e., their 
reactive energy densities vanish. Hence the reactive energy of 
their sum, the standing wave, is due entirely to the interference 
between the two traveling waves. That is, the invariants in 
(3) consist only of the cross-terms. Furthermore, since every 
globally sourceless field is a Fourier superposition of null plane 
waves with uo = ck, it follows that the reactive energy of every 
globally sourceless EM field is due entirely to self -interference. 
This gives a partial intuitve explanation of EM rest energy, 
as seen most clearly in the standing wave example. However, 
the rest energy of fields with sources need not be entirely due 
to self-interference since their Fourier synthesis also requires 
plane waves with uj ^ ck, which are not null. (Such plane 
waves represent 'virtual photons,' which have positive mass.) 



II. The energy flow velocity v(r, t) 

The correspondence between the rest energy Eq of a relativis- 
tic point particle and the reactive energy density 7Z(r,t) of an 
EM field in vacuum can be extended to include the velocity 
of the point particle, 



E ' 



whose field-theoretic version is 



v(r,t) 



cS(r,t) 



(7) 



(8) 



U(r,t) ' 

Poynting's theorem d t U + cV • S = — J • E then becomes 

dtU + V-(vU) = -J-E, (9) 



which shows that U behaves like the density of a compressible 
fluid with source —J E, flowing at velocity v(r, t). Note that 

v(r,t) = \v(r,t)\ = c <^> K(r,t) = 0. (10) 

Thus, while the field (E,B) propagates at c, its energy 
generally flows at v < c almost everywhere. 

For the standing plane wave of Example 2, (4) shows that 
\v\ < c as expected, and 

v = <=> cos 2 (kz — ujt) — cos 2 (kz + ujt) 
<^> kz + out = ±(kz — out) + nir. 



v = <^> z 



tn (11) 



Hence v has fixed nodes in both space and time: 

fin nir 
— — = z n or t = — i 
2k 2uj 

where n is any integer. Since v changes sign at z n and t n , the 

energy is totally reflected at these nodes. 

The energy oscillates back and forth between the nodal planes 
z = z n , and v(z,t) = z ■ v oscillates between d=c at any z. 

The conflict between the moving nodes (6), where v = ±c, 
and the stationary nodes (11), where v = 0, is resolved by 
noting that v(z, t) is undefined when U = and S = 0, so 

cos(kz — out) — cos(kz + uot) = 0. 

This gives cos kz cos uot = and sin kz sin ujt = 0, hence 

z = z n and t = t rn with m + n odd. 

These planes are the intersections of the traveling and sta- 
tionary nodes. Intuitively, the reason why v(z,t) is undefined 
at these events is that perfect reflection there requires it to 
change instantaneously between the values ±c. At all other 
values of z, v(z,t) still oscillates between d=c but does so in 
a continuous manner; see Figure 6 in [1]. 

III. A HISTORICAL NOTE 

The fact that the energy of an EM field in vacuum generally 
flows at speeds less than c was noted almost a century ago 
by Bateman [3, page 6]. To the best of my knowledge, 
this important insight has remained undeveloped and largely 
unappreciated. 2 I believe this phenomenon, and its relation to 
reactive energy as detailed in [1], are fundamental features 
of electromagnetic fields which ought to be studied both 
theoretically and experimentally. 
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2 For a time-harmonic field of frequency uj, the energy transport velocity is 
commonly defined as Voj(x) = cSuj(x)/Uuj(x), where Suj(x) and U u (x) 
are the time averages of S(r, t) and U(r,t) over one period 2ir/uo. In general, 
l v a;(2c)| < c almost everywhere. However, time-averaging is lossy and the 
ratio of two averages is not the average of the ratio. Hence v w is not a 
time average of the exact, instantaneous energy flow velocity v(r,t). I thank 
Professor Andrea Alu for pointing this out. 



